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^ ■ Abstract 
00 . 

t^T* ■ The anisotropic Bianchi I cosmological model coupled with perfect fluid is quan- 

tized in the minisuperspace. The perfect fluid is described by using the Schutz 
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like equation is solved and a wave packet is constructed. The norm of the resulting 
wave function comes out to be time dependent, indicating the loss of unitarity in 



• formalism which allows to attribute dynamical degrees of freedom to matter. A 

, Schrodinger-type equation is obtained where the matter variables play the role of 

I time. However, the signature of the kinetic term is hyperbolic. This Schrodinger- 

O 

$—1 . 

this model. The loss of unitarity is due to the fact that the effective Hamiltonian 
is hermitian but not self-adjoint. The expectation value and the bohmian trajecto- 
ries are evaluated leading to different cosmological scenarios, what is a consequence 
of the absence of a unitary quantum structure. The consistency of this quantum 
d I model is discussed as well as the generality of the absence of unitarity in anisotropic 

quantum models. 

PACS number(s): 04.20.Cv., 04.20. Me 

1 Introduction 

One of the main hopes regarding quantum cosmology is the possibility to obtain the initial 
conditions that will determine the ulterior evolution of the Universe when its classical 
regime is reached [1, 2]. The task of obtaining a quantum cosmological scenario faces 
many difficulties, one of them being the absence of a natural time variable, since the 
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general relativity action leads to a constrained system which is invariant under time 
reparametrization. There are many attempts to recover the notion of time in quantum 
cosmology [3]. For example, a time coordinate may be identified with the space volume, 
which is a growing function in an expanding universe. But, all these attempts have 
revealed of limited applications until now, and the problem of time in quantum cosmology 
remains an unsolved puzzle. 

Another suggestion to incorporate a time variable in quantum cosmology is through 
matter fields. This proposal has been extensivelly discussed in reference [4]. It has been 
shown that a typical ordinary quantum mechanical structure can be built up: a Hilbert 
space, with an inner product, as well as sets of physical observables, may be identified. 
All analysis performed in reference [4] was made keeping the functional character of the 
Wheeler- De Witt equation. A simpler consideration in the same sense was made in refer- 
ences [5, 6, 7, 8, 9, 10, 11], where the matter fiuids were introduced with the aid of the 
Schutz's variables [12, 13], the Wheeler-DeWitt equation being written in the minisuper- 
space. The employement of the Schutz's variables permits again to identify the matter 
fields with time, since the associated momentum appears linearly in the Lagrangian; the 
restriction to minisuperspace has the advantage of allowing an explicit intregration of the 
resulting Schrondinger-like equation. 

Although of phenomenological nature, these quantum cosmological perfect fiuid mod- 
els in the minisuperspace are a very good laboratory in order to verify the consistency 
of constructing quantum cosmological models where the time variable is identified with 
the matter fields. Since a Schrondinger-like equation is obtained, all the machinery of 
ordinary quantum mechanics can be employed. In the references [5, 6, 7, 8, 9, 10, 11], this 
has been done for isotropic universes in many different situations, connected mainly to the 
nature of the matter content. In order to keep the effective Hamiltonian hermitian, the 
inner product in the Hilbert space has acquired an aditional factor. Boundary conditions 
on the wave functions were imposed, assuring the hermiticity (and, at the same time, the 
self-adjointness) of the effective Hamiltonian operator. Wave packets were constructed 
from which specific predictions were obtained by computing the expectation value of the 
observables (in this case, the scale factor) or by evaluating the bohmian trajectories. 
Since the modulus of the wave function integrated in all space is time independent, both 
results agree. The main conclusion of those works is that the quantum model predicts 
a singularity-free universe which exhibits a bounce approaching the classical behaviour 
asymptotically. 

In the present work we will attempt to extend this analysis to anisotropic cosmological 
models, specifically to Bianchi I models. The minisuperspace approach will be used, as 
well as the description of the matter fields through the Schutz's formafism. The initial 
aim is to verify if quantum effects may suppress the anisotropics in the same way they 
have suppressed the initial singularity in the isotropic case. In doing this analysis, an 
unexpected feature of anisotropic quantum cosmological model appears: the norm of the 
wave function comes out to be time dependent. Hence, the quantum model is non-unitary. 
This leads immediately to the question if this is a real quantum system. With a non- 
unitary theory a probabilistic interpretation can not be implemented, since the norm of the 
wave function is not a conserved quantity. At the same time, the ontological interpretation 
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of Bohm-de Broglie becomes doubtfull, since bohmian trajectories are not conserved also. 
In both cases, some kind of "creation" must be admitted in order to have some quantum 
interpretation of the rcsuhs. If an anisotropic quantum cosmological perfect fluid model 
is a legitimate quantizable system, some fundamental changes in quantum mechanical 
interpretation must be implemented. 

The reason for this loss of unitarity when a simple extension from isotropic to anisotropic 
universes is made, is due, in our opinion, to the fact that the kinetic term of the effective 
Hamiltonian is not positive definite and, at the same time, due to the measure in the 
original gravitational action. The hyperbolic signature of the kinetic term implies that 
there is a whole line in the phase space where the momenta are not zero, and yet the 
energy is zero. Hence, the wave function need not to be zero at infinity along this line, 
and this implies that the probability current is non-zero at least at some points at infinity. 
Hence the time derivative of the norm of the wave function is not zero anymore. It is 
important to notice that the effective Hamiltonian is hermitian. The loss of unitarity 
comes from the fact that this effective Hamiltonian is not self-adjoint and it does not 
admit any self-adjoint extension. 

This is not an artifact of the construction of the wave packet, but a general feature of 
anisotropic quantum models which lead to an hyperbolic signature to the kinetic term of 
the Hamiltonian function with an unusual measure. This suspiction is supported by the 
fact that when the kinetic term is made elliptic by force, the norm of the wave function 
becomes time independent and normal quantum framework is restablished. But, the 
hyperbolicity of the Hamiltonian is not the only reason for the lack of unitarity: the 
specific measure in the action, due to the determinant of the metric, is also essential in 
the sense that its suppression by force restore also the unitarity independently of the 
signature of the kinetic term. Hence, the loss of unitarity is a direct consequence of a 
gravitational system (which leads to a non conventional measure in the inner product) 
which exhibits anisotropies (which leads to a hyperbolic kinetic term) . 

As it will be verified, the loss of unitarity leads to a fundamental discrepancy between 
the many-worlds intepretation (based on the Copenhagen interpretation) [14] and the 
Bohm-de Broglie interpretation [15, 16] of quantum mechanics. In both cases the universe 
exhibit a bounce. But the computation of the expectation value for the metric functions 
reveals a universe always isotropic; on the other hand, the bohmian trajectories reveal a 
universe where anisotropies are present near the bounce disappearing asymptotically. It 
is important to stress that such quantum Bianchi I model has been extensivelly studied 
in the literature [17, 18, 19, 20, 21]. But, in all previous work no matter field has been 
used, and hence no explicit time coordinate has been identified. For this reason, in our 
opinion, the loss of unitarity has not been remarked before. 

This article is organized as follows. In the next section, wc construct the Wheeler- 
DeWitt equation for the anisotropic perfect fluid model and we determine the wave func- 
tion by using the separation of variables method. In section 3, a wave packet is constructed 
and its norm is shown to be time dependent. The reason for this unexpected result is 
discussed. In section 4, the expectation value of the scale factor and the bohmian trajec- 
tories are obtained. The discrepancy between them are settled out. In section 5 we discuss 
the results and present our conclusions. In the appendix we show how the equivalence 
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between many-worlds and dBB interpretation disappears due to the absence of unitarity 
in a quantum model. 



2 Wheeler-DeWitt equation for an anisotropic per- 
fect fluid model 

Our starting point is the action of gravity coupled to a perfect fluid in the Schutz's 
formalism: 

A= f d^x^f^R + 2 I d^xVhhabK"'' + [ d^x^f^p (1) 

Jm JdM JM 

where K'^^ is the extrinsic curvature, and hab is the induced metric over the three- 
dimensional spatial hypersurface, which is the boundary dM of the four dimensional 
manifold M; the factor IGvrG is made equal to one. The first two terms were first ob- 
tained in reference [2]; the last term of (1) represents the matter contribution to the total 
action in the Schutz's formalism for perfect fluids, p being the pressure, which is linked 
to the energy density by the equation of state p — ap. In the Schutz's formalism [12, 13], 
the four-velocity is expressed in terms of five potentials e, /9, and S: 

u.^-{e,, + cP,, + es,,) (2) 

A* 

where /i is the specific enthalpy. The variable S is the specific entropy, while the potentials 
C and P are connected with rotation and are absent for FRW's type models. The variables 
e and 6 have no clear physical meaning. The four velocity is subject to the condition 

= 1 . (3) 

The metric describing a Bianchi I anisotropic model is given by 

ds' = N'^df- (x{tfdx'' + Y{tfdy^ + Z{tfdz''\ . (4) 



In this expression, N{t) is the lapse function. Using the constraints for the fluid, and after 
some thermodynamical considerations, the final reduced action, where surface terms were 
discarded, takes the form 



A= J dt - (^XYZ + XZY + YZX^ 



(5) 



At this point, is more suitable to redefine the metric coefficients as 
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Using these new variables, the action may be simphfied further, leading to the gravita- 
tional Lagrangian density 

Prom this expression, we can evaluate the conjugate momenta: 

g3/3o g3/3o p'i(^o 

Po = -12^/3o , P+ = 12^/5+ , p_ = 12— /3_ . (8) 

The matter sector may be recast in a more suitable form through the canonical trans- 
formations 

T^pse-^p;^''+^^ , PT^p':^'e^ , e^e-{a + l)^ , = . (9) 

Pe 

The final expression for the total Hamiltonian is 

i/=7Ve-3/3o|_l(p2_^^_^2j^g3(l-a)/3o^^| _ (10) 

The fundamental aspect of the Hamiltonian (10) to be remarked is the hyperbolic signa- 
ture of the kinetic term. 

The lapse function N plays the role of a Lagrangian multipher in (10). It leads to the 
constraint 

H^O . (11) 

The quantization procedure consists in considering the Hamiltonian as an operator which 
is applied on a wave function 

= (12) 
taking at the same time the momenta as operators (we use natural units where ^ = 1): 

Since the momentum associated to the matter degrees of freedom appears linearly in the 
Hamiltonian, we can identify it with a time coordinate 

= . (14) 

Due to the canonical transformations employed before, this new time is related to the 
cosmic time tby rft = e3"/3oc?T. In this way, we end up with the Wheeler- De Witt equation, 
in the minisuperspace, for an anisotropic Universe filled with a perfect fluid: 

92 92 92 



^^_24ie=^a-«)A^ . (15) 



The wave function must obey the following boundary conditions: 



^'Ift^ioo = «^|ft^±oo , (16) 

with K e (— cxd.cxd], /3j denoting the dynamical variables. These boundary conditions 
are established by requiring that the Hamiltonian be hermitian. For k = and oo, the 
boundary conditions are ^''|ft^±oo = and ^'|/3j^±oo = 0, respectively. As it will be 
seen later, in spite of being hermitian, the effective Hamiltonian is not self-adjoint and 
does not admit any self-adjoint extension. This will lead to the loss of unitarity. A more 
detailed discussion on the sclf-adjoint properties of the operators in quantum cosmology 
with perfect fluid, in a situation very close to the present one, is given in references [7, 11]. 
A rigorous mathematical discussion is given in reference [22]. 

Now, our goal is to solve (15) and to construct the corresponding wave packet. To do 
so, we use the separation of variable's method. First, we write the wave function as 

^{Po, P+, T) = cj>{f3o, (3+, (5_)e-'^^ , (17) 

leading to the equation 

The function is then written as 

0(/3o,/3+,/3-)=To(/3o)T+(/3+)T_(/3_) , (19) 

leading to the equation 

^ + 24^e3(i-«)A - ^ - ^ = (20) 
To T+ T_ ^ ^ 

where we have simplified in an obvious way the notation for the partial derivatives. The 
solutions for the functions T± are 

T± = CiC^'^i^i , (21) 

where C± are constants and k± are the separation parameters. These separation param- 
eters must be real otherwise the wave function is not normalizable. 
The equation determining the behaviour of Tq takes then the form, 

To"+(^24£;e'(^-")* + (A;^ + A;!))To = , (22) 

the primes meaning derivatives with respect to Pq. It is easily to see that the parameter 
E must be positive. The previous equation can be solved through the redefinitions 

= 6^^° , y = a' , r = ^{l-a) . (23) 
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after what (22) takes the form of a Bessel's equation: 



To + — + 

y 



To f2AE 1 



Tn = 



(24) 



where A;^ = fci + k'^_ and the dots are derivatives with respect to y. The solution is 



\ r 



V r 



with u — ik/r, Ci,2 being integration constants. 

The final expression for the wave function is then 



= e 



-iET 



(25) 



(26) 



where Ci,2 are combinations of the preceding integration constants. 



3 The wave packet: Loss of unitarity 

We want now to construct a superposition of the solutions (26), generating a regular 
wave packet. In principle, this can be achieved by considering the integration constants 
as gaussian functions of the parameters k± and E. The general case constitutes a hard 
problem from the technical point of view. We may consider, for simplicity, the final 
wave function independent of one of the variables which amounts to fix one the 
corresponding parameters k^ or k^ equal to zero. From here on wc will consider A;_ = 0. 
Notice that the final results would be the same if we had imposed A:+ = and k_ ^ 0. 
Hence, even if the anisotropic models are not analyzed in all their generality, a large class 
of them is covered in what follows. 

Fixing A;_ = 0, the wave packet is given by 

* = /e^'=+/^-{(7iJ.(^a^)+(72J_.(^a^)}e-'^^dM£^ . (27) 

In principle, in the expression for v it appears the modulus of A';+ while in the first 
exponential in (27) we have — oo < A;+ < +oo. We will consider a superposition of both 
Bessel's functions in such a way that the expression for the wave packet may be written 
as 

* = A{k+,q)e'^+'^+J,[qa'^e-''''^dk+dq , (28) 

with q = ^^^^ and 

A{k+,q) = e-^'''iq''^'e~^'^' . (29) 
In this case, the integrals can be explicitly calculated, leading to the wave packet 

a^'- {p+ + C{a,B)f- 




4:B 47 



(30) 
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where 

B = X + isT , C(a,B) = lna-—^ — .\n2B , s = -^il^l^ . (31) 

3(1 — a) 32 

Notice that the wave packet given by (30) is square integrable, and it vanishes in the 
extremes of the interval of vahdity of the variables a — and except along the hne 
/3o = — /9+ where it takes a constant value, being consequently regular as it is physically 
required. The wave packet (30) is indeed a solution of the equation (15), as it can be 
explicitly verified, and it obeys the boundary conditions fixed before. If we discard the 
terms corresponding to the variable /5+ (connected with k+), the wave packet for the 
isotropic case [10] is reobtained. 

The main point to be remarked now is that the norm of (30) is time dependent. Using 
the definition a — and integrating in and a we obtain 

r r a'-'-^*^dadp^^-f^^^F{T) . (32) 
Jo J-oo 3(1 — a) \ 



where 



and 



F(T) = exp(g) , (33) 



C{a,B) = CR + iCi , 

CR = \na-—^\nAB*B , Ci = ^ . ~^ . arctan(^) . (34) 
3(1 — a) 3(1 — a) A 

The norm of the wave function is time dependent. Hence, the quantum model is not 
unitary. 

The absence of unitarity may be understood by inspecting again the wave packet (30). 
In fact, this wave packet goes to zero at infinity, excepted along the line = — /?_|_, where it 
takes a constant value at infinity. This does not spoil the regularity of the wave packet; in 
particular, it remains finite when integrated in all space and specific boundary conditions 
are obeyed. But, this leads, at the same time, to an anomaly in the infinity boundary. 
The reason for this anomaly may be understood by analysing again the Schrodingcr-like 
equation (15). Notice that, after decomposing it into stationary states, the energy E is 
zero along the whole line f3o = — Along this line, the wave function need not to vanish. 

It would be expected that a hermitian Hamiltonian operator should always lead to 
a unitary quantum system, since the Hamiltonian operator is responsible for the time 
evolution of the quantum states. The problem here relies on the fact that, in spite of 
being hermitian, the Hamiltonian effective operator 

H,ff^e-'^'-'^^[d',-dl-dl] (35) 

is not self- adjoint. This means that W = H but the domain of is not the same as the 
domain of and the conservation of the norm becomes senseless [23]. 
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In order to verify if an operator is self-adjoint or not, we must compute the so-called 
deficiency indices n± which are the dimensions of the linear independent square integrable 
solutions of the indicial equation 

H(j) = ±i(p . (36) 
Using the effective Hamiltonian (35), the solutions of the indicial equations are ^ 

(f)+ = CiJ„{y) + C2J-uiy) , (37) 
cf>. = c^K^iy) + cj,{y) , (38) 



where u — ik/r and y = ^Ji/r'^a^, a and r having the same definitions as before. It 
is easy to see that J±v{y) and Iv{y) are not square integrable solutions while K„{y) is. 
Hence, n+ = and n_ = 1 and, as explained in [23], the effective Hamiltonian operator 
is not self-adjoint and does not admit any self-adjoint extension. Notice that changing 
arbitrarily the signature in (35) or suppressing the unusual measure, the deficiency indices 
become n+ = n_ = and the effective Hamiltonian becomes self-adjoint. 

It is important to remark also that, due to hyperbolic character of the Hamiltonian, 
the energy E may take negative values. Hence, we may expect that this system is unstable 
since the energy is not bound below. However, it is possible to consider a kind of " Dirac 
sea" hypothesis, with all negative energy state filled, and to take effectively into account 
only positive energy states. 

To support the idea that this anomaly is due to the hyperbolic signature in the kinetic 
term in (15), together with the unusual measure, let us change it to an elliptic signature 
by force. In doing so, the main change in the wave functions (55) is that the order of the 
Bessel function becomes real: u — \k\. We keep only the Bessel function of positive order 
because it does not diverge as a — > 0. In evaluating the norm of the wave function, we 
consider the same superposition factor as in (29) and we first integrate on the parameter 
q, obtaining. 



-fe+(c{a,B)+ip+y 



dk . (39) 



Now we writte 

N 



/-\-oo roo 
/ a^-^'''^*^dadp+ , (40) 
-oo ^0 

and we integrate first in f3^ and then in a. The unusual measure in the integrals is due 
to the requirement that the reduced Hamiltonian in (15) must be hermitian [10, 11]. The 
final result is 

N^^'l 2^/T(l + fc/r)e-'^*^' , (41) 
Jo 

where T{x) is the gamma function, and is a new constant. The norm of the wave 
function is finite and, more important, time independent. The same occurs if instead the 
measure is suppressed. 

^It must be remarked that in finding the solutions, we supposed the system to be independent of the 
variable /?_ and we perform a plane wave expansion in the variable In this sense, we considered just 
a one-dimensional system depending on the variable (3q with a parameter k. However, the conclusions 
do not depend on these considerations and we could consider at least a two-dimensional system by, for 
example, performing a gaussian superposition in the parameter k. 
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4 The scenario for the Universe 



The fact that the quantum cosmological perfect fluid model leads to a non-unitary quan- 
tum system implies in principle that no usual quantum interpretation can be be applied 
to it, unless we allow creation of universes. In what follows we will adopt the point of 
view that this is a legitimate quantum system which ask for a convenient framework in- 
terpretation. Hence, we will try to extract previsions for the evolution of such a universe 
using the many-worlds and dBB interpretations scheme. Of course, these interpretations 
scheme must be enlarged in order to incorporate non-unitary quantum system. It is not 
sure that this can be done consistently. However, this very important conceptual problem 
is outside the purpose of the present work. Here, in particular wc will show that the 
many-worlds and Bohm-de Broglie interpretations leads to different results. This is due 
to the lack of unitarity as it is explained in the appendix. 

Before to do this, let us just recall the classical solutions for the Bianchi I cosmological 
model with a barotropic perfect fluid described by p = ap. For the time parametrization 
dt — a^°'dT, t being the cosmic time, the functions X, Y, and Z admit the solution 

X{T) = e*+^++^^-=Xo(r + c)'^(r-c)'^ , (42) 

Y{T) = e'^°+'^+-^''-=yo(r + c)'^(T-c)'^ , (43) 

Z{T) = e^°-2'^+ = Zo(r + c)'^(r-c)'^ , (44) 

where c is constant, and si, S2 and S3 are parameters such that 

si + S2 + S3 = , + + = 6 . (45) 

Notice that there is an initial singularity, near which the Universe is very anisotropic, 
becoming isotropic asymptotically. 

Let us return now to the computation of the quantum scenario through the use of the 
many-worlds and ontological interpretations of quantum mechanics. 



4.1 Expectation values of the dynamical variables 

Given the wave function the expectation value of a variable /3i is obtained in the usual 
way: 

^ e^^^-<^)^o^*^dPodP^ ■ ^ ' 

For /3j = /5o in (46) we find for the numerator: 
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where we have noted 



n 



exp(—ti) In ti (iti ~ —0.577 , u 



A 



2B*B 



3(l-a) 



Hence, 



This result leads to 



3(1 - a) 



In 



2\B\ 



1 + 



A2 



3(l-c«) 



(48) 



(49) 



(50) 



where gq is a constant. This is the same result as in the isotropic case [10]. Consequently, 
the space volume evolves as in the corresponding isotropic case. 

The anisotropics are represented by the function whose expectation value will be 
computed in what follows. We will evaluate now the numerator of (46) with /3j = 
Integrating in and expressing (3q in terms of a as before, we find: 



/+00 r 
-oo J—oo 



2B*B 

The integrals Ii and I2 take the form. 



^"exp 



r HAB*Bl IF{T) 



In a da 



^"exp 



A 



3(l-a) 



1 



9 1 - a 



A 



n + In 



2B*B 
A 



/2 = 



25*5 



25*5 



3(1 -a) A 



da (52) 



(53) 



We find finally 



,(l_„).--v .-, . (54) 

The expectation value of /3+ does not depend on time. Consequently, the predicted result 
for the evolution of the Universe in this case is the same as in the isotropic case: there is 
no anisotropy during all the evolution of the Universe. A similar computation shows that 
< > is also time independent. The cosmological scenario is really isotropic. 



4.2 Computation of the bohmian trajectories 

The result found in the last section indicates no trace of the anisotropics existing in the 
classical model in the corresponding quantum analysis. We will evaluate the bohmian 
trajectories which determine the behaviour of a quantum system in the ontological inter- 
pretation of quantum mechanics. 

In the ontological interpretation of quantum mechanics, the wave function is written 

as 

* = E exp(ie) , (55) 
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where S is connected with the amphtude of the wave function, and to its phase. 
When (55) is inserted into the Schrodinger's equation, the real and imaginary parts of the 
resulting expression leads to the conservation of probability and to a Hamilton- Jacobi's 
equation supplemented by a term which is identified as the quantum potential, which 
leads to the quantum effects distinguishing the quantum trajectories from the classical 
ones. 

In this formulation of quantum mechanics, the trajectories (which are real trajectories) 
corresponding to a dynamical variable q with a conjugate momentum pg are given by 



Pi 



dq 



(56) 



The ontological formulation of quantum mechanics leads to a natural identification of a 
time coordinate, what is very important for quantum cosmology where in general there 
is no explicit time coordinate. 

Let us consider the wave function (30). Putting in the form (55), the phase reads. 



e(/3o, T) = - arctan(^ ^ ^ , ^^^^ 



(57) 



where all quantities are defined as before. The conjugate momenta associated to the 
dynamical variables (3o and (3+ read 



(58) 



where we have explicitly used the time parametrization such that the lapse function is 
given hy N — a^". The bohmian trajectories are then given by the expressions 



3(1 — a) , 1. „ a 



3(l-a) _ 



'AB*B 27 ' 

"27 ' 

dots representing derivatives with respect to T. Combining (59,60), we find 

sT 



-12a^-'^"d = 3(1 - a) 



4:B*B 

This last equation leads after integration to the expression 



+ 



(59) 
(60) 

(61) 



a e 



3(l-t«) 



(62) 



Reinserting the relation in the equations (59,60) we can obtain the following solutions to 
a and 



/ — 1 \ 3(l-c) 

V24sA7y 



3(l-a) 



/34 



1 



3^^_^^ln{arctan^(f^)+E 



arctan^ 



sT 



+E 



A 

- 24sA7 



3{l-a) 



3(l-a) 



lnL> 



(63) 
(64) 
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where E and D are integration constants. Remember that s < 0. 

In opposition to the expressions obtained for the expectation values of /Sq (which is 
connected to a) and f3^ in the preceding subsection, the bohmian trajectories predict an 
anisotropic Universe. Until this point, this strange discrepancy is not so catrastrophic: 
in order the bohmian trajectories coincide with the results for the expectation value for 
some quantity, the integration constants that appear in the former must be averaged over 
an initial distribution given by the modulus of the wave function at T = 0. At T = 0, we 
have 



/3+(r = 0) = In 



-24sA7D3(i-a) 



E 



Hence, 



7^, = ***|r=o 



TT ^ 48s7 18(l-c«)^7 



In^ 



3(l-a) 



.3(l-a) 



(65) 
(66) 

(67) 



For Pq and the average on the initial conditions leads to the integral expressions 

_ fOO fOO 

Po{T) = / / e'^'-^^f'on, Po{T) dPi, dPi , (68) 
Jo Jo 



PAT) 



Jo 



g3(l-a)/3j^^ /3+(T) dpi, d(3\ 



(69) 



where f3^ and f3\_ denote the initial values of the metric functions. In the isotropic case 
[10] the expression corresponding to the above ones leads to a perfect agreement between 
many-worlds and dBB interpretations. These expressions can be recast in the following 
form: 

1 ^1,2 . y ^^^^^^ (1j.^ ,^.)2 

:i - q;)2 -24s7A io Jo 



Po{T) 



9(] 



ln|(A 



2 ^ ^2y2N 
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(70) 



(71) 



where x = 4~3{i-a)/} and y = E. The variables x and y were restricted to positive values 
in order to assure that the metric functions are real. Even if the integrals (70,71) seem to 
admit no simple closed expressions, it is evident that they are time dependent. 



5 Conclusions 



It is generally expected that quantum effects in the very early universe may furnish the 
set of initial conditions which will determine the subsequent evolution of the Universe 
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when its classical phase is reached. By initial conditions we mean here the isotropy and 
homogeneity. Moreover, it is also expected that those quantum effects may lead to the 
avoidance of the initial singularity, one of the major problems of the standard cosmological 
model. In this work we have tried to analyse the possibility that quantum effects can 
suppress initial anisotropies. Specifically, we have studied a Bianchi I model with a perfect 
fluid, with an isotropic pressure, employing the Schutz's description for perfect fluids. This 
problem has for us two main interests: flrst, it adds more degrees of freedom with respect 
to the isotropic model, since now we have four independent variables instead of just two; 
second, it permits to verify if anisotropies in the early Universe disappear in the quantum 
model, as it happens with the initial singularity for the corresponding isotropic one. The 
employement of Schutz's formalism for the description of the perfect fluid present in the 
model allows us to identify quite naturally a time coordinate associated to the matter 
degrees of freedom, since the canonical momentum corresponding to the matter variables 
appears linearly in the Hamiltonian. Hence, the Wheeler-DeWitt equation can be reduced 
to a Schrordinger-like equation in terms of three dynamical variables related to the metric 
function, /3o, /?+ and /?_. In order to treat the problem analytically, we have restricted the 
problem to the special case that the wave function is independent of one of the variables, 
namely 13 -. 

The resulting Schrodinger-like equation has a hyperbolic signature in the kinetic term. 
This means that a state of zero energy is possible along an inflnite line where po = P+ 
or, equivalently, (5q — This leads to an anomaly in the boundary at inflnity which, 
however, does not spoil the regularity of the wave function. But as consequence, the 
resulting quantum system is not unitary anymore. The reason for this loss of unitarity is 
the absence of a self-adjoint extension for the hermitian effective Hamiltonian (35). The 
lack of self-adjointness in this model is due both to the hyperbolic signature (which is 
a consequence of treating an anisotropic quantum cosmological model) and due to the 
presence of a non-trivial measure in the original Hamiltonian (10) (which is consequence 
of treating a gravitational system), as a detailed inspection of the computation of the 
deficiency indices reveals. The loss of unitarity leads to a discrepancy between the many- 
worlds and Bohm-de Broglie interpretations of quantum mechanics. Notice that if wc go 
back to the isotropic model, where there is no hyperbolicity anymore, but yet an unusual 
measure, the unitarity is restored. It is important to remark that, in our opinion, the loss 
of unitarity pointed out here in anisotropic models has no relation with the description 
for the matter fluid, since the hiperbolicity of the kinetic term and the unusual measure 
appear already in the pure gravitational sector. 

In order to verify this explicitly, we have solved the Wheeler-DeWitt equation in the 
minisuperspace. A wave packet was constructed being regular in the sense that it is 
square integrable. Using this wave packet, wc have determined the behaviour of the 
metric functions using flrst the many-worlds interpretation of quantum mechanics, which 
implies to compute the expectation value of those functions. We found that there is 
no trace of anisotropies at any moment: the expectation value of the function is 
constant while the expectation value of (5q has essentially the same expression as in the 
isotropic version of this problem. All the features of this model are the same as in the 
isotropic case. Later, we have determined the behaviour of metric functions employing the 
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ontological interpretation of quantum mechanics, determining the bohmian trajectories. 
In this case the function is no longer a constant, and an initial anisotropic Universe 
is predicted. Asymptotically, it becomes isotropic like in the classical case. This result is 
maintained even after the averaging on the initial conditions. As it is well known [15, 16] 
the bohmian trajectories should lead to the same results that are obtained computing the 
expectation values after averaging on the initial conditions. This equivalence does not 
occur for the anisotropic Bianchi I cosmological model because the equivalence between 
both interpretations is valid only for unitary system. 

Let us now precise in another way that the loss of unitarity is due to the anomaly 
at the boundary at infinity generated by the hyperbolic signature of the kinetic term. 
Writting the wave function as in (55), we obtain the expression 



which should express the conservation of probability. Integrating in /?o and /3+ does not 
lead to the vanishing of the integrated left hand side, since the norm of the wave function 
is time dependent. The term in the right-hand side can be converted to a surface integral 
at infinity. The currents are zero at infinity excepted along the line /3o = and hence 
the right hand side also does not vanish. 

The question resulting from the analysis made above is if it is possible to consider 
seriously such quantum cosmological perfect fluid model. In principle, the answer, in our 
opinion, must be positive, since this system is a natural extension of the corresponding 
isotropic case, where everything is well definite. But, in order to take seriously such 
model, the usual interpretation scheme of quantum mechanics must be enlarged in order 
to take into account creation/annihilation of probabilities and bohmian trajectories. In 
what concerns cosmology this perhaps may be done since these creation/annihilation refer 
to disconnected universes. But, this remains an open issue, and what we may assert at 
the moment is that the loss of unitarity is a general feature of quantum system where the 
signature of the kinetic term is hyperbolic as it happens with anisotropic cosmological 
models. Moreover, since the predictions for the evolution of the universe using the many- 
worlds or the ontological interpretations do not coincide, one of them must be more 
suitable for implementation of this enlarged interpretation of quantum mechanics. Notice, 
finally, that the loss of unitarity appears already when quantum fields arc quantized in 
space-times with closed timelike curves [24]. But this represents another context since 
the space-time itself remains classical. 
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for many enlightfuU discussions and CNPq (Brazil) for partial financial support. 
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A Conditions for the equivalence between many- worlds 
and dBB interpretations 

Let us for simplicity consider a one dimensional quantum mechanical system: 
Prom this expression we obtain 

where we have used (55) and we have defined TZ — ^*^. The expectation value and the 
averaged bohmian trajectories are given by 

< X >c= '^7S ,xb= I TZox{t) dxo (75) 

J_^l-cax J-oo 

where C and B stand for "Copenhagen" and "Bohm", respectively. The subscript in 
the second integral indicates that the quantities must be evaluated at t = 0. Also in 
this second integral the function x{t) is obtained by integrating the bohmian trajectories 
x{t) = dQ/dx. If the norm of the wave function is made equal to one in t = 0, those 
quantities are identical at this moment. 

To show that they are the same for all time, we must just show that their derivatives 
are the same for any value of t. Taking the derivative and using the expression for the 
current in the first integral and the expression for the bohmian trajectories in the second 
one, we obtain 

d /-"^ „ 90 , d „ de 



TZ — dx , —xb = TZq — dxQ , (76) 

-oo ox dt J-oo ox 

where we have assumed the norm of the wave function equal to one. Integrating the 
expression for the current, we have that the derivative of the norm of the wave function 
is zero. Hence, it has the same value at any time, from which 

/oo roo 
TZtdxt = / TZodxo . (77) 
-oo J —oo 

This implies that, after a changing of variable in the second integral, 

7^o^ = 7^* . (78) 

dxt 

Inserting this in (76) wc find that the derivative of expectation value and of the averaged 
bohmian trajectories are the same for any value oft. Since, they were equal for t = 0, both 
quantities are identical for any value of t. Notice that it was essential to have a constant 
norm of the wave function in order to obtain this result. In the Bianchi I quantum model 
studied in this paper the analysis follows the same lines, but the norm of the wave function 
is time-dependent and the equivalence exhibited here is no longer valid. 
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